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Abstract 

We propose a non-standard subsampling procedure to make 
formal statistical inference about the business cycle, one of 
the most important unobserved feature characterising fluc- 
Qh' tuations of economic growth. We show that some character- 

. istics of business cycle can be modelled in a non-parametric 

way by discrete spectrum of the Almost Periodically Cor- 
■ related (APC) time series. On the basis of estimated char- 

acteristics of this spectrum business cycle is extracted by 
Ph I filtering. As an illustration we characterise the man prop- 

' erties of business cycles in industrial production index for 

Polish economy. 
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dure. 
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^ 1 Introduction 

O " 
oo . 

Tj- , Seminal works, that originated interest in empirical modelling of business cycles in 

I macroeconomy, clearly indicated theirs inhomogeneity for both, spatial a nd time domain . 



O I In particular, introductory remarks of W.C Mitchell in Business Annals. iMitchelll ( 19261 ). 

contains the following suggestion: No two recurrences in all the array seem precisely alike. 
Business cycles differ in their duration as wholes and in the quickness and the uniformity 
with which they sweep from one country to another. When identifying business cycle R.E. 
Lucas proposed its own conception, which, as he states in his 1977 paper, identifies the 
business cycle with movements about trend in gross national income. These movements 
are typically irregular in period and in amplitude. Regula rities are onl y observed in the 



co-movements among different aggregative time series; see iLucasI (119771 ). It is clear, that 
Mitchell initially suggested different time pattern of business cycles for different economies. 
However, it is obvious, that from the dynamic point of view, as Lucas states, business 
cycle exhibit irregular and nonperiodic character. 
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For develop ed economies som e stylis ed facts about busines s cycle s are known in the lit- 
erature; see iKing and Rebelol fll999h or IStock and Watsonl fll999h . But, we see the lack 
of precise and well established methods of formal statistical modelling of those empirical 
properties. It prompts new studies resulting many different approaches and frameworks 
of business cycle extraction: see for example exhaustive review presented Diebold and 
Rudenbush fll996h . When the lack of the theory of statistical inference seems to be a per- 
sistent state, the consensus about empirical properties of business cycles is based, either 
on an ad-hoc reasoning, or on the empirical results, that are possible to confirm using a 
group of methods, built on the basis of relative different frameworks. However, the ex- 
traction of the business cycle component from observed time series is still a controversial 
issue. In particular, since there is ongoing interest in many approaches to separate growth 
component from the cyclical component, and because there is no consensus on how to 
detrend the data, the business cycle stylised facts are sensitive to the adopted procedure. 
Hence, this has become not only a controversial issue in the business cycle theory itself, 
but also a subject of criticism by competing empirical approaches, as well. 
The main purpose of this paper is to present a novel approach to formal business cy- 
cle estimation. We propose a non-standard subsampling procedure, in order to make 
formal statistical inference about the properties of the business cycle. We show, that 
business cycle can be modelled by parameters of discrete spectra of the Almost Periodi- 
cally Correlated (APC) stochastic process. The AFC class is a g eneralisatio n of Fe riodi- 
cally Correlated (FC) class of stochastic processes, introduced by iGladyshevI (Il96ll ). The 
vast literature confirmed substanti al empirical importance a n d flexibility of FC c lass i n 
many time series app l icatio n s, see: iFarzen and Fagand (119791). lOsborn and SmithI (119891 ). 
Franses and BoswiikI (Il996n. iFra nses iBoUerslev and Ghvselsl (I19961). Franses and 

Ooms (ll997lV iBurridge and Taylor (2001). iFranses and DiikI (120051). According to Hurd 



and Miamee (120071 ). the periodically correlated stochastic processes are nonstationary, but 
non-constant unconditional expectation of the process exhibit periodic, and hence regular, 
evolution in time domain. The generalisation presented in this paper assumes that the 
mean of the nonstationary time series can be described by almost periodic function, i.e. 
the function, that belongs to the topological closure of periodic class of functions. 
From the definition, AFC stochastic processes may describe irregular character of uncon- 
ditional means for nonstationary time series. Assuming, that detrended time series follows 
AFC, we relax assumption of stationarity of cyclical factor, very commonly imposed in 
filtering approaches. Nonstationarity of the cycle component of the series, together with 
possible irregularities in time pattern of the unconditional mean, makes our approach rel- 
atively flexible and general. Consequently, incorporating the AFC factor into the model 
of observed discrete time series should result in much more accurate approach to business 
cycle extraction than those proposed so far. 

In the empirical part of the paper we analyse the cyclical behaviour of production sector in 
Foland with the use of the model with AFC stochastic component. We characterise busi- 
ness cycle on the basis of industrial production index and also on some subsector indices. 
We discuss the empirical results and reflect them to the previous analyses conducted for 
the Folish economy. 
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2 Basic notation and definitions 



Formally, the second-order and real-valued time series {Xt : t G Z} is called periodi- 
cally correlated if both the mean function n{t) = E{Xt) and the autocovariance function 
B(t, t) = cov{Xt, Xt+r) are periodic at t for every r G Z, with the same period T. In order 
to introduce the class of almost periodically correlated time series we need the definitio n 
of almost periodic function. We recall the following definition from ICorduneanul ( 119891 ): 



Definition 2.1. A real-valued function f{t) : Z — y M of an integer variable is called 
almost periodic (AP in short), if for any e > there exists an integer > 0, such that 
among any consecutive integers, there is an integer with the property 

sup|/(t + p,)-/(t)| <e. 

A second-order real-valued time series {Xt : t G Z} is called almost periodically 
correlated if both the mean function fi{t) = E{Xt) and the autocovariance function 
B{t,T) = cov {X t, Xt-i-r) are almost periodic function of an integer variable, for every 
r G Z. It is easy to see that any periodic function is almost periodic. Therefore, the class 
of APC time series is wider than the class of PC time seri es. Duri n g last five decades 
the APC cl ass was broadlv applied in telecommunication (Gardner! (19861). Napolitano 
and SpoonerlioOll)), climatology iBloomfield et al.l (|l994l)') and ma ny ot hers fields. For 



exhau stive review of possible applications see iGardner et al.l (l2006l ) and ISerpedin et al. 
( 120051 )). Empirical importance of such a class of nonstationary time series prompted new 
studies concerning properties and estimation methods. 

In APC case the mean function and the autocovariance function B[t, r) for any r G Z 
has the Fourier representation of the form: 



B{t, 



T ~ 



AeAr 



a(A, r)e 



iXt 



(1) 



where the Fourier coefficients m{ilj) and a(A, r) are given by: 



1 

m(^) = lim — \ jji{t)e~ 
t=i 



(A,r)= lim i^i?(j,r)e- 



iXj 



(2) 



see 



Hurdl ( ll99lh . iHurd and Miamed (l2007h l According to ICorduneaniJ (jl989l ) sets 
= {tp G [0,27r) : mx{ip) ^ 0} and A^ = {A G [0, 27r) : a(A,r) 7^ 0} are count- 
able. Hence, the set A = {J^^^Aj. is also countable. If the time series is PC, then 
representations ([T]) become equations and the sets and A are contained in the set 
{2kTT/T : A; = 0,1,...,T- 1}. 

In the problem of business cycles extraction the vast econometric literature exploit ap- 
proaches based on the assumption of zero mean imposed on the distribution of stochastic 
factor describing business fluctuations. Moreover, this stochastic factor is usually mod- 
elled under stationarity assumption, leading to the framework that utilises parameters 
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of continuous spectrum. The econometric approach presented in this paper relaxes sta- 
tionary assumption, and consequently a more general dynamic model of observed time 
series is subject to empirical analysis. We model business cycles in a non-parametric way, 
taking into account discrete spectra of observed time series. It means that we characterise 
business cycles by non-zero frequencies ip E and by corresponding Fourier coefficients 
m (ib). The definition and properties of discrete spectra in simple representation see Priest- 
lev (1l98lh . or in PC case in lHurd and Miamed (120071 ). 



Notice that any ipQ e (0, 27r) corresponds to the length of the cycle 27i/ipQ. Hence the 
following testing problem: 

Ho: V^o ^ ^ /ON 

Hi: V^oe^, ^"^^ 

enables to test the statistical significance of the cycle with appropriate length. According 
to the definition of the set and Fourier coefficients m(-) our testing problem is equivalent 
to the following: 

Hq: |m(V'o)|=0 
Hi: |m(V^o)|^0. 

We consider formulation H] in details. Since we are interested in business cycle estimation 
we restrict frequency ipo such that corresponding length of the cycle is not shorter then 
1.5 years. This formally means, that in further analysis for monthly data we assume that 
V^o e (0,0.35). 

In this paper by amplitude, which corresponds to frequency G \I' fl (0,0.35), we define 
the distance between maximum and minimum value of the function h{t) = 2Iie[m{ip)e^'^*]. 
The problem stated above requires statistical theory of detecting significant frequencies 
in the set and corresponding Fourier coefficients m(-). In the next sections we present 
some results given sampling model generated by APC assumption. 



3 Estimation problem 



By {Xc,^+i, Xc„+2, • • • , ^c„+dn} "we denote a sample from APC time series {Xt : t G Z}, 
where {rfrijneN is any sequence of integers tending to infinity with n and {c„}„gN is any 
sequence of integers. For any ip G [0,27r), estimator m'^''^"{ip) (m^''^(?/^) for short) of the 



parameter m{'ijj), in representation 
the form: 



[T]), based on sample {Xc„+i, X, 



Cn-\-dn 



-iipj 



c„+2,---,^c„+d„} takes 



(5) 



The standardi sed version of ([5]) has asymptotic normal distribution with zero mean; see 



LenartI (j2011bl ). Theorem 2.1. Additionally, the variance of this distribution is a function 
of values of the generahsed spectral density, calculated at arguments dependent on tjj. 
For de f inition and basic proper ties of generalized spectral density in AP C case s ee iHurdl 
(ll989h . lDehav and Hurdl (Il994l )l Since the standard theory, presented in iLenartI (1201 lal ). 
provides methods of estimation of generalised spectral density in APC case only under 
the zero-mean assumption or under the assumption that the set \E' is known and finite. 
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estimation of generalised spectral density in our case is not possible so far. Therefore, 
in the paper we exploit subsampling methodology, to construct asymptotically consistent 
test related to (jlj). In this approach the asymptotic variance estimation is not of particular 
in terest. Simila r ly, su bsampling methodology was also used for PC case in time domain 
in lLenart et al.l (120081 ) and for APC case in frequency domain in iLenartI (l2011al ). 
The problem of frequency es timation ijjo ca n be solved on the basis of a more generalised 
approach than presented by IWalkerl (jl97ll ). Given assumption that there exists interval 
/^p, such that n \E' = {ipo}, it is possible to formulate the natural estimator of the 
unknown frequ ency ipo of the form ipn = arg maXa;g/^^{A/n|m„(x)|}. As it was shown in 
LenartI ( 1201 Ibl ). Theorem 3.1, under some regularity conditions we have: 





P ^ 


m(V'o) 


'4'n 







(6) 



4 Subsampling procedure and consistency 

In this section we describe the main id ea of the subsamplin g methodology, according to 
the approach presented and developed by iPolitis et al.l (119991 ). We use the same notation. 
Initially we assume that the time series {Xf : t G Z} is governed by unknown probability 
distribution P, that belongs to a certain class of probability measures V. Denote by 
{Xi, X2, . . . , Xn} a sample from the time series {Xt : t G Z}. Our goal is to approximate 
the distribution of 

i;„(ft„-^(P)), (7) 

where 9n = On{Xi,X2, . . . ,Xn) is an estimator of 9{P), the parameter of interests, and 
Vn is appropriate normalising sequence. Let b{n) {b for short) be any sequence of integer 
numbers tending to infinity with n, such that b < n and b/n 0. 

One of the main assumption in subsampling methodology is that there exists asymptotic 
distribution of ([7]). We denote this distribution by J(P), with J{x, P) as a corresponding 
cumulative distribution functions at point a; G M. Following the idea of iPolitis et al.l ( 119991 ) 
the distribution of ([7]) can be approximated by its subsampling version of the form: 



Ln,b{x) 



1 



n-b+l 



n 



— ^ l{Vh{9n,b,t - On) < X}, 



(8) 



t=l 



where is the indicator function of the set B and = ^b(Xt,Xt+i, . . . ,Xt+5_i) 

as an estimator of the unknown parameter 9{P) obtained on the basis of the sample 
. . . with if: as a starting point and & as a siz e of su bsample. Under 

suitable regularity conditions it is known that (see IPolitis et al.l (119991 ). Theorem 4.2.1, 
page 103): 

(i) if X is a continuity point of J(-, P), then Ln h{x) J{x, P), 

(ii) if J{-,P) is continuous, then sup3..g]K \Ln,b{x) — J(x,P)| — ^ 0, 
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(iii) If J(-, P) is continuous at point c(l — a), then 



P (vn{9n - e{P)) < c„,fe(l - a)) ^ 1 - a, (9) 

where for any a G (0, 1), we define 

c„,6(l - a) = inf{x : Ln,b{,x) > 1 - a}, 

c(l — a) = inf{x : J(x, P) > 1 — a}. 

The imphcation (iii) is crucial to construct a subsamphng confidence interval for the pa- 
rameter 9{P). 

We are interested in estimation of the absolute value of coefficients of the Fourier repre- 
sentation of the mean of APC process. Namely, we take 9{P) = \m{ip)\. S ubsamp l ing pro - 
cedure, with 9n,b,t = '^n = y/n, is consistent; see for details LenartI ( 2011b ). 

Theorem 2.3. Consequently, the confidence intervals for the parameter 6{P) = \m{il))\, 
obtained by subsampling procedure, are asymptotically consistent. 

Now let take any ipQ G (0,7r]. The test (jlj) with test statistics Unii'ip}) = A/n|^n(V')| and 
subsampling critical value are asymptotically consistent. In our paper we prove some mod- 
ification of this result (see Theorem 18.21 in Appendix). We use test statistics n„({?/'}) = 
y/n\fn{4')\, that can be interpreted as a value of test statistics Iln{{i'}) = y/n\mn{4')\ 
based on the sample {Xi — X„, X2 — X„, . . . , X„ — X„}, where X„ is the sample mean 
for the path {Xi,X2, . . . ,Xn}. The critical value c^^i^ ~ calculated according to 
the formula that utilises subsampling procedure: 



ci?(l-«) = inf{x:LX^(x) >!-«}, 

n-b+l 

n — b + 



n-b+l 

n + i 



where 



n 

" J=c„+1 



e^(^) = - (^.-X„)e-^^-. 



and fnijp) = "T^'" ('?/'). 



5 Statistical model of cyclical fluctuations 

In this section we present the statistical framework of extraction the cyclical com- 
ponent when the one-dimensional time series describing changes in economic activity is 
observed. In Section 15.11 we present basic assumptions concerning the model, while in 
section we describe in details the algorithm of formal statistical extraction of business 
cycle component. 
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5.1 Model structure and assumptions 



Let consider a real- valued time series, denoted by {Pt : t G Z}. At the beginning of 
this section we assume that the unconditional expectation for the process {Pt : t E 1^} 
exists for any t E'L. 

An interesting case, that is of particular interest in econometrics is the class of integrated 
stochastic processes, denoted by I{d) for integration of order d G N. If we are interested 
in analysis of I{d) processes in our framework, some additional assumptions should be 
imposed top assume the existence of unconditional moments. In the case when {Pt : t G Z} 
is J(l) process it is sufficient to assume additionally that there exists to E Z such that 
E{Ptg) < oo. Hence we obtain in this case, that E{Pt) < oo for any t E 'L. More 
generally, if {Pt : t G Z} follows I{d) process, then it is sufficient to assume that there 
exists to £ ^! such that E{Ptg+k) < cxd for /c = 0, 1, . . . , ci — 1, to assure moment existence. 
Consequently, we formally exclude in our analysis processes with pure integration, but 
some restricted cases, representing strict nonstationarity with finite unconditional mean, 
may be modelled. 

For further analysis we assume that the mean function ^p{t) = E{Pt) is defined by the 
sum of deterministic function f(t,(3), parameterized by /3 G MP, and almost periodic 
function g{t), with the Fourier expansion of the form: 



git) = J2 mp^e^^'. 

For convenience, we rewrite g{t) in equivalent representation: 

g{t) = ap{il)) cos{il)t) + hp{il))sm.{il)t). 

i/;e<I'pn[0,7r] 

This automatically implies, that: 

/ip(t) = /(t, /3) + g{t) = fit, mpii^y'^'. 



(10) 



Equation ( ITTj) leads to a more general approach to modell ing business fluctuations, tha n 
those presen te d in the literature so far; s ee for example: iBeveridee and NelsonI (Il98ll ). 
ClarkI fll987h . iHarvev and Jaegej fll993h . iHamiltonl fll989h . iKrolzieJ ( 1l997h . The main 
advantage of our approach is, that it only relies on the specification of the first moment 
of the time series {Pt : t G Z}, making model assumptions much weaker. To illustrate the 
importance of our assumptions and generalisation we present an example below. 

Example 5.1. Let {Pt : t G Z} be a time series such that Pt = Pt-i + et, where 
E{Pq) = b and {ej : t G Z} is APC time series with expectation function /Xe(-) such 
that fie{t) = a + g{t) — g{t — 1), where : Z — )■ M is a function of the form: 



git) 



m{ip)e 



iipt 
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a G M and card(\I/) < oo. Notice that for any t > 1 we have 



= Po + ei + £2 + • • • + e*. 

Therefore 

t 

E{Pt) = b + J2 E{ej) = b + at- g{0) + g{t) = f{t, /3) + g{t), 
i=i 

where f{t, (3) = f3o + f3it, (3o = b — g{0), (3i = a. This means that time series {Pt : t E Z} 
can be represented as ([TT]) . If g(t) = 0, then He(t) = ci, and time series {Pt : t G Z} can 
be interpreted as /(I) process with drift and assumption E{Pq) = b. 

The function f{t, (5) can be interpreted as a trend component, modelled in this paper 
by the polynomial. The function g{t) contains summarised information about seasonal 
fluctuations, business fluctuations and long-term cyclical fluctuations. From the Fourier 
representation of g{t) we split the whole set \Efp of non-zero frequencies into the mutually 
exclusive sets, that are related to those three periodic attributes of time series dynamics. 
Initially, we interpret long-term cyclical fluctuations as those with the length more than 
8 years, since the frequency uj is related to the length of cycle that equals 271 /uj units. In 
order to distinguish cyclical fluctuations from seasonal fluctuations we assume formally, 
that in the representation ffTOl) the set = {ip '■ m{ip) 7^ 0} C [0,27r), is unknown. For 
the set "^p, let consider the following decomposition: 

= ^p_i U ^P,2 U ^'p^s. (12) 

We assume that "^pi fl (0,0.35) = ^Pp^i, and consequently the set \E'p^i represent all 
frequencies with corresponding length of the cycle greater than 17 months. Therefore 
the set \E'p,i contains frequencies that can describe business fluctuations. The set \E'p,2 
contains only seasonal frequencies, namely \E'p^2 C {2fc7r/12 : k = 0,1, . . . , 11} while \E'p_3 
contains all remained frequencies. In the following section we concentrate our attention 
only to parameter identiflcation and estimation in the set ^E'p,!. 

5.2 Cycle identification and estimation 

Our approach aims at identiflcation and estimation of cyclical fluctuations. In order 
to remove trend component and to weaken seasonal effects, observed time series is subject 
to some preliminary transformations, . Hence, we formulate the algorithm of frequency 
identiflcation that consists of three basic steps. The flrst step enables to remove sea- 
sonal component, the second step detects the trend component, while in the third step, 
parameter identiflcation and estimation is provided, 
ep 1 - removing the seasonal component. Seasonality appear in most monthly eco- 
nomic time series. More formally, we allow (it is assumed), that \E'pn{2A;7r/12 : k = 
1, 2, ... , 11} 7^ 0. Since, the estimation of the frequencies and corresponding Fourier 
coefficients, that represent seasonal frequencies is not of particular importance in 
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our paper, we use centered moving average filter 2 x 12MA (see: iMakridakis et al. 



fll998h . iBrockwell and David f l2002h ) to remove seasonal pattern. We show below 
that this filter does not change the elements of the set "^p^i, what is crucial for future 
estimation procedure. Denote by {Yf : t G Z} time series obtained by application of 
the centered moving average filter. It means, that = L2xi2{B)Pt, where 

L2MB) = {B-^ + 25-^ + . . . + 2B-^ + 2 + 25 + . . . + 25^ + B^)/2A, 

and B^Pt = Pt^k for any t and k. Note that the expectation of the time series 
{Yt : t G Z} exists. On the basis of the Theorem 18.11 and elementary calculations 

we get 

/xy (t) = E{Yt) =Po + Pit + ... + l3pt^+y2 ^Y{^)e'^\ (13) 

/(i,/3) 

where n {2fc7r/12 : A; = 1, 2, . . . , 11} = 0, ^y = \ {2A;7r/12 : A; = 1, 2, . . . , 11} 
and / is a function. Fourier coefficients mp{ip) and mY{ip) are related according to 
the formula: 

my(V') = L2xi2(e-*'^)mp(7/>). (14) 

Notice that /(t, (3) is also a polynomial of order p. In particular, for s G {p — l,p} 
we have Ps = /3s- Consequently, given model with p = or p = 1 (i.e. constant or 
linear trend) we have, that fit, (3) = fit, (3). In case p = 2 functions / and / have 
different values, but fit, (3) — fit, (3) is constant over time. Additionally, filtering 
the series with centered moving average operator, we obtain, that: 

^y n (0,0.35) = ^p,i 

and 

n^P,2 = 0, 

which means that the set \l/y still contains the same elements as \E'p,i and dos not 
contain seasonal frequencies from the set \E'p,2. 
Step 2 - removing the trend component. The case when p = is trivial. Let consider 
the case p = 1. Application of difference operator -Z^i(-B) = (1 — i?) for the time 
series {Yt : t G Z} results with time series {Xj : t G Z}: 

X, = L,iB)Yt = Ft - Yt_, = (P,+6 - P,_5 + Pt+5 - Pt-6)/24. 

The expectation of the time series {X^ : t G Z} exists and is described by almost 
periodic function of the form: 

/xx(t)=/3i+ Yl ^xi^)e'^\ (15) 

where 

C {0} U ^p \ {2A:7r/12 : A; = 1, 2, ... , 11}, (16) 
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and 

n (0,0.35) = ^p,i, 
which follows from the Theorem 18. 1[ Additionally, we have: 



(17) 



:i8) 



and the Assumption 18.11 holds. 

In the general case, when p G M we use natural operator Lp{B) = (1 — By. The 
resulting time series {Xt : t G Z} can be represented by the following transformation 
of Yt. 

Xt = il- BfYt, 



and hence, the expectation of Xt takes the form: 

E{Xt) = ^lx{t) 



(19) 



where, according to the Theorem 8.1: 



Lp{e 



-iil>\ 



mp{i)). 



(20) 



By estimation \Lp{e ^^)L2x12{g ^^)\ > 0, which is true for any p G N and G 
(0, 0.35), we have: 

n (0,0.35) = *P,i. (21) 

Therefore the problem of parameter identification and estimation in the set \l/p_i 
reduce to the problem of parameter identification and estimation in the set H 
(0,0.35). 

Step 3 - parameter identification and estimation. The formula fl2Tl) is crucial in the 
algorithm of parameter identification and estimation in the set \l'p,i. Initially, in 
Step 3 we formulate the additional assumption that the autocovariace function of 
time series {Xt : t G Z} exists and it is almost periodic function. Notice, that 
the weaker assumption concerning periodic structure of autocovariance function 



Parzen and Paeano 


(1979). 


Osborn and Smith 


(1989). 
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1996). 


FransesI fl996h.lFranses and Ooms 


(1997). Franses and Dijkl (2005h. We use statistics 



Iln{{ip}) = ^/n\fn{ip)\ and corresponding critical value Cn,fe(0.99%) for the series 
generated from the previous steps of the algorithm. The test statistics Iln{{ip}) 
can be interpret as a value of test statistics !!„({?/'}) = ^/n\'mn{'^p)\ based on the 
sample {Xi — X„, X2 — X„, . . . , X„ — X„}, where X„ is the sample mean for the 
path {Xi, X2, . . . , Xn}. The critical value is ca lculated acco r ding to the formula 



that utilises subsampling procedure presented in iPolitis et al.l ( 119991 ): 



a) = mf{x:Glf,\x)>l-a}, 
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where 



1 



n-b+l 



n-b+1 



t=l 



We fix 6 = 2.'b^fn and we calculate test statistics and corresponding critical value 
for i\} from the discrete set of frequencies on the interval (0, 0.35). If the value of test 
statistics is greater than the critical value on some subinterval J C (0, 0.35) we take 
this subinterval as the interval containing some elements of the set p^\. Next, we 
estimate the frequency connected with subinterval / using ^ . By plug in technique 
we estimate amplitude related to each identified frequency in almost periodic part 
of the mean function of the process \Pt : t G Z}. 
In our algorithm step 3 is fundamental in procedure of extraction business cycle component 
from the observed time series. Its main advantage is, that the frequencies describing the 
cyclical dynamics of economic activity are subject to formal statistical inference. This 
clearly distinguishes our approach from many other procedures presented in the literature, 
where the lack of statistical uncertainty in the procedure is ve ry comrnon an d for ces ad-hoc 
appro ach; see for example polemics concerning detrending in lCanoval ( 119981 ) and iBurnside 
fll998h . 

However, it is very important, that the procedure yields only statistically significant 
frequencies, and extraction of the business cycle is subject to additional filtering. In the 
empiric al part of the paper we use the Hodrick-Prescott fHPl filter fsee Hodrick and 



Presco t t (119971)'), with smoot h ness parameter A. According to iGomezI ( 119991 ). [Gomez 



( l200lh . iMaravall and del Rid (120011 ) parameter A can be described as the argument of 



frequency V'o 



A 



(22) 



4(l-cos(V^o))'' 

where 2tt /ipQ can be interpret as a length of the cycle. Hence, on the basis of our procedure, 
it is possible to choose appropriate parameter A of the HP filter, restricting spectrum only 
to significant parameters in the set p^i. Alternatively it is possible to apply any filter in 
cycle extraction. We choose the simplest HP filter for illustrative purposes. 



6 Empirical illustration 

In this part of the paper we analyse cyclical behaviour of production sector in Poland. 
In particular we apply our model and three step procedure in order to characterise business 
cycles in industrial production index and some subsector indices. 

Figure [T] (a) presents time series of industrial production indexQ in Poland from January 
1995 to December 2009 (2005 year = 100%). This index contains: mining and quarrying; 
manufacturing; electricity, gas, steam and air conditioning supply. In the first step we 
applied centered moving average filter 2 x 12MA to eliminate strong seasonal effects. The 
results of filtering is plotted on the Figure [T] (b). It is clear, that centered moving average 
filter removes the seasonal effects and also business fluctuations are clearly observable (see 

^Source: Eurostat. 
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Figure [11(b)). 

According to our algorithm, presented in previous sections, we present in Figure [2] (a) 
first differences of tlie centered moving average filter applied for industrial production 
index. We see some evidence about the existence of cyclical behavior in time series under 
consideration. The amplitude of cycle does not seem to be constant over time. Also, the 
amplitude is smaller in period 1995-2001, while after year 2001 is characterized by greater 
variability. Therefore we use logarithm transformation for industrial production index 
to stabilize the amplitude. Figure [2] (b) presents the first difference of centered moving 
average filter applied for logarithm of industrial production index. It is easy to see that 
the amplitude is more constant over time then before logarithm transformation. 
Figure [3] presents plots of the values of the test statistics Iln{{ip}) = \/n\fn{ip)\ with 
corresponding critical value c„^b(0.99%). The test statistics exceeds the critical value in 
three subsets on the interval (0, 0.35). Hence, taking care only of significant values of test 
statistic on the Figure |3] and in zoom on Figure HI we assume that: 

^P,in(0,0.35) = {V'i,V^2,V'3}. (23) 

The values of estimated frequencies from the set "^pi were calculated according to (E]). 
These values with corresponding length of the cycle can be found in Table 1. 
Estimated amplitude of the cycle with corresponding length 8.5 and 3.4 years equal 0.13 
and 0.07 respectively. This second amplitude dominates the estimated value of amplitude 
with corresponding length 2 years. We see, that 8.5-year length of the cycle received 
data support. However we can not formally interpret such fluctuations as a long-term 
growth trend or business fluctuations. We should rather look at this fluctuations as a 
mixture of both long-term growth trend and business fluctuations. Consequently and 
unquestionable, the dataset support fluctuations with corresponding length 3.4 years as 
a basic characteristic of business cycle in industrial production in Poland. To confirm 
this statement we extract cyclical fluctuations from industrial production index (filtered 
by centered moving average filter 2 x 12MA ) with the use of HP filter condition to the 
values of parameter A fixed for A = 5 500, A = 12 000, A = 32 000, A = 55 000. The results 
are plotted on Figure |5l By restricting parameter A to values stated above, according 
to the formula |22l we extract fluctuations with the length not greater than 4.5, 5.5, 7 
and 8 years respectively. Since our goal was to extract only business fluctuations without 
significant influence of long-term growth trend, we restrict filtering only to fluctuations 
with corresponding length shorter than 8 years. 

Analysing plots presented on Figure |5] it is possible to confirm the presence of cycles in 
industrial production in Poland with estimated length in the interval 3-4 years (during the 
period 1995 - 2009). In Table [2] we determined the periods of recessions and expansions 
in industrial production. We interpret turning points as margins of this periods. In most 
cases the recession is shorter than expansion. Consequently, our analysis confirm results 
discussed in the literature, that business cycle in industrial production for Poland display 
asymmetric behaviour. Also, the business cycle t roughs are rather sharper than peaks, 



which is also typical for business cycles; see iHicksl ( ll950l ). iMilas et al.l (120061 ) 



In the next step we provide a more detailed analysis based on a formal identification of 
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business cycles in sectors and subsectors of industrial production in Poland. We use the 
same statistical tools as for the total industrial production index. We considered all cat- 
egories identified for industrial production. The set of all modelled indices are presented 
in Table |3l 

Figure [H] presents plots of logarithms of all considered indices. Repeating the procedure, 
initially applied for the total index, we use centered moving average filter 2 x 12MA to 
remove seasonal effect from the data sets (see Figure [7]). First differences are presented 
on the Figure [HI It is clear, that majority of indices exemplify cyclical pattern, just like 
in the case of index of total production, but with rather differential amplitudes and length. 
To identify frequencies in the unknown set \E'p,i we applied again the test statistic Iln{{ip}) = 
V^\''^n{'ip)\ and corresponding critical value clfjjia)- The results are presented on the Fig- 
ure [9], where we plotted estimated lengths of the cycles together with appropriate esti- 
mated amplitude. In different sectors and subsectors of industrial production the data 
provide evidence in favour of cycles with length in the interval 1.5-3 years. However, 
those cycles are characterized by much shorter estimated amplitude than the cycles with 
length in the interval 3-4 years. The cycles with estimated length in the interval 3-4 years 
were supported in predominant set of subindices. Only in the case of manufacture, food 
products and beverages (C10_C11), manufacture of basic pharmaceutical products and 
pharmaceutical preparations (C21) and electricity, gas, steam and air conditioning supply 

(D)- 

We see, that the largest estimated amplitude characterize cycles with the corresponding 
length of more than 4 years. But only in a few cases the cycles with length 5-8 years 
were supported. It can be seen clearly on the Figure [TTl where the comparison of all 
identified cycles for all 32 indices is presented. In spite of the fact, that observed time 
series were subject to filtering with the use of centered moving average filter 2 x 12MA , 
all investigated subindices provide data support in favour of the existence of cycles with 
length not greater than two years. However, as seen on Figure [TTl those short cycles are 
characterized by amplitudes with values located relatively close to zero, as compared with 
longer significant cycles. This makes such a short term periodic pattern not extremely im- 
portant in describing cyclical behaviour of modelled time series. Additionally, all indices 
support cycles of length 3-4 years, with relatively greater value of corresponding ampli- 
tudes as compared to characteristics of short term fluctuations. Also, except manufacture 
of wearing apparel (14-th index) we see no data support for cycles with length between 4 
and 7 years. Consequently, for all considered subindices, the set of statistically significant 
cycles is clearly divided in two separate parts. The first set is constituted by short term 
cycles with small amplitudes together with middle term fluctuations, attributed in most 
cases by stronger amplitudes. The second set consists of frequencies, describing long term 
cycle, namely with length not less than 7 years. Just like in case of the total production 
index, we tend not to interpret those long term fluctuations as important characteristic of 
business cycle for Polish economy. According to our results, just like for the total index, 
all considered subindices are characterized by existence of the long term trend. 
Using HP filter we extract business cycles from all industrial production indexes. Similar 
as for industrial production index - total we fix the parameter A as A = 5 500, A = 12 000, 
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A = 32 000, A = 55 000 (see Figure [T2|) . The reasons why we chose those values of A 
parameter are the following. Firstly, we fix the same parameter to compare results with 
those obtained for industrial production index. Secondly, the length of the cycle that is 
greater than 8 years is not clearly constant over different subindices and therefore we can 
not interpret those fluctuations as business fluctuations. Finally, we can notice that in 
the interval from 4 to 8 years there are only a few significant lengths of cycles and this 
should give rise to extract similar shape of business fluctuations for different values of 
parameter A. Almost all extracted fluctuations reveal presence of cycles with length in 
the interval 3-4 years. Summing up, the cycle with length in the range 3-4 years is typical 
and prevalent for cyclical fluctuations in industrial production in Poland. 



7 Concluding remarks 



In this paper a novel approach in business fluctuations analysis for one dimensional 
economic processes is proposed. Using theory of almost periodically correlated time se- 
ries and subsampling procedure we consider a formal approach to estimate the length of 
business cycles. The main advantage of our approach is, that the business cycle charac- 
teristics are treated in formal way, and are subject to statistical inference. This clearly 
distinguishes presented framework from many filtering-based approaches, broadly consid- 
ered in empirical applications. We model business fluctuations by parameters of discrete 
spectra of time series, under assumption that amplitude of this fluctuations is constant 
over time. Taking in consideration estimated length of the cycles we extract business 
fluctuations by HP filter for parameter of smoothness chosen on the basis of formal pro- 
cedure. 

The main conclusion presented in empirical illustration is that, during period 1995-2009, 
we confirm (using statistical tools) the presence of 3-4 years length of business cycle in 
industrial production index in Poland. This result was obtained either on the basis of 
the total index and also analysing subindices. This result co nfirms analyses conducte d 



so far on the bas i s of Polish macroeconomic time series: see iGra dzewic z et al. fl2010f ). 



Adamowicz et al.l (120081 ) , ISkrzypczyhskil (120081 ) and ISkrzypczyhskii (.2006i ) . 
All indices and subindices supported significance of short term and middle term fluc- 
tuations, attaching relatively small amplitudes for periodicity with length less than 2 
years. Additionally, in all time series we detected existence of longer term cycle (7-8 
years), interpreted in this paper as a trend or a mixture of both trend and business cycle 
fluctuations. 



8 Appendix 

Assumption 8.1. Let {Xt : t E "Z} be APC time series such that for any x G [0,27r) 
there exists a constant B{x) (which dependents only on x), such that we have estimation 



E 

V'e*\{a;} 



m{ilj)cosec 



X 



< B(x) < oo, 



(24) 
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Theorem 8.1. Let {Xf : t G N} he a time series for which the expectation function exists 
and it is almost periodic function of the form fix{t) = E{Xt) = X] We 

assume that for the set \I/ and corresponding Fourier coefficients m(-) the Assumption \8.1\ 
holds. Let L{B) = Yl'j=-p^jB'' linear filter, where p,q > 0, {aj}j=_p is a sequence 
of real numbers, and B^Xt = Xt-j for any j G Z. Then 

where \l/y = andniyi'ip) = L{e~'^^)mx{4') ■ Additionally, assumption \8. 1\ holds for the 
set \l/y and corresponding coefficients myi-)- 

Proof. Notice that 



\j=pi 



P2 



P2 



j=pi Ve*x 



j=pi 

ip&^x j=Pl ip&^x 



(25) 



P2 



By estimation |my(?/')| < \mx{4')\ S kil 'we conclude that condition 1.1 from iLenart 



0=P1 



f 2011bl ) holds for the set \E'y and corresponding Fourier coefficients 



my( 



□ 



Theorem 8.2. Take any ip G (0,27r). Let the assumptions of Theorem 2.2 in iLenart 
Eollh ) hold. Then 

(i) L^^}{x) A J^^\x), for any x G M, 

(ii) sup,,^|LX>(a;)-jW(a:)|-^0, 

(iii) subsampling confidence intervals for the parameter \m{ilj)\ are asymptotically con- 
sistent, which means that 



p(v^i\f^m-\mm)<cSa~c.) 



1 — a, 



(26) 



where b = b{n) — )■ oo and b/n — t- 0. 



Proof of t he TheoremjS.^ I this proof we use the same steps as in Theorem 4.2.1, 
page 103 in IPolitis et al.l ( 1999 ). Let g = n — 6+1, r„ = ^/n and 
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Notice, that 

= \ilnn[\K-''\n-\rnm]+n[{\mm-\rn^^^^^ < 

^ t=l 

We need the following lemma. 

Lemma 8.1. For any real x and e > we have estimation 

Unix - e)l{E4 < L;^^/(x)1{E4 < Unix + e), (27) 
where = {n max |(|m(^/')| - \fni^)\) + i\fi~^'\^)\ - |m^r^'''(^)l)l < e}- 

l<t<q 

Proof. Let consider two cases: 

1° l{En} = 0, inequality ([27]) holds 
2° l{En} = 1, then 

r, max |(|m(^)| - |f„(V^)|) + (le^'^WI - < e, 

l<t<g 

which means that for any I <t < q 

vi-''\^) ■■= niiimin - \rnim + iK-''\n - \K-''\m] e 

Using next inequality x — rj^j^^^'^i'ip) > x — e, which is true for any 1 < t < g we get 

Hn[\mi-'\^)\ - \min] + n[i\mi^p)\ - |f„(v^)i) + nc'^n - \K-''\m] < ^} 

= l{r4K-i'^(^)| - |m(^)|] <x-v'-''\^)} > l{n[\mi-'\n ' I^WI] <x-e}. 

(28) 

Analogically, using inequality x — Vj^n^'^iip) < x + e we get 

nn[\K-'\n - \rnim + niiMn - \rnim + iK-'\n - \<'''\m] < ^} 

= l{n[\rn}-''\n - |m(^)|] <x- r,l^'\ij)} < l{n[\m'~'\n - |m(V^)|] < x + e}. 

(29) 

Summing inequality ( 125]) . fl29|) dla t = l,2,...,g we get f l27|) . This completes the 
proof of lemma. 

□ 

In next step we show that P(-E'n) 0. Using inequality — |2;2|| < \zi — -22! (which 
is true for any complex numbers zi,Z2) and inequality | ^j=p Cje*-?^] < Cp|cosec(x/2)|, 
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(which is true for any x ^ modulo 27r and real numbers Cp > Cp+i > • • • > Cg ) we have 



max |r/^ ^'^(^)| < max n\\m{ilj) \ - |f„(V^)|| + max n\\fl^ - ^ 

< n\\m{i,)\ - |r„(^)|| + max rfelf^-i'^(^) W I 

l<t<q 

< n\m{ij) - f„(^)| + ^|X„||cosec(^/2)| (30) 

< n\m{ilj) - rhn{ip)\ + Tb\mn{ip) - r„(^)| + ^|X„||cosec(z/'/2)| 



< rb\m{ip) - mn{ip) \ + ^|X„||cosec('?/'/2)| + -^|X„||cosec(^/'/2)| 

By convergence Tb\m{ilj) — m„('?/^)| A and y|X„| A we get 

max|r/^i'''(^)| AO, 



which means that P{En) 1. Using next Slutskys Lemma and Theorem 2.1 in iLenart 
( 2011b ) we have \/n{\r^~'^'^{il))\ — |m(-?/')|) -4 To finish the pro of it is sufficient to 



follows next steps in Theorem 4.2.1, page 103 in iPolitis et al.l (119991 ). therefore we omit 



them. □ 
References 

E. Adamowicz, S. Dudek, D. Pachucki, and E. Walczyk. Synchronizacja cyklu koniunktu- 
ralnego polskiej gospodarki z krajami strefy euro w kontekcie struktury tych gospodarek. 
Instytut Rozwoju Gospodarczego, Szkoa Gwna Handlowa, Warszawa, 2008. 

S. Beveridge and C.R. Nelson. A new approach to decomposition of economic time series 
into permanent and transitory components with particular attention to measurement 
of the business cycle. Journal of Monetary Economics, 7(2):151-174, 1981. 

P. Bloomfield, H. Hurd, and R. Lund. Periodic correlation in stratospheric ozone time 
series. J. Time Ser. Anal, 15(2):127-150, 1994. 

T. Bollerslev and E. Ghysels. Periodic autoregressive conditional heteroscedasticity. Jour- 
nal of Business and Economic Statistics, 14(2):139-152, 1996. 

P.J. Brockwell and R.A. Davis. Introduction to Time Series and Forecasting. Springer- 
Verlag, New York, 2002. 

C. Burnside. Detrending and business cycle facts: A coment. Journal of Monetary 
Economics, 41:513-532, 1998. 

P. Burridge and A.M. Taylor. On regression-based tests for seasonal unit roots in the 
presence of periodic heteroscedasticity. Journal of Econometrics, 104:91117, 2001. 

F. Canova. Detrending and business cycle facts. Journal of Monetary Economics, 41: 
475-512, 1998. 

P.K. Clark. The cyclical component of U.S. economic activity. Quarterly Journal of 

Economics, 102:797-814, 1987. 
C. Corduneanu. Almost Periodic Functions. Chelsea, New York, 1989. 



17 



D. Dehay and H. Hurd. Representation and estimation for periodically and almost pe- 
riodically correlated random processes. In: W.A. Gardner (Ed.), Cyclostationarity in 
Communications and Signal Processing, IEEE Press, pages 295-329, 1994. 

F.X. Diebold and G. Rudenbush. Measuring business cycles: A modern perspective. 

Review of Economic and Statistics, 78:67-77, 1996. 
P.H. Pranses. Stochastic Trends in Economic Time Series. Oxford University Press, New 

York, 1996. 

P.H. Franses and H.P. Boswijk. Temporal aggregation in a periodically integrated autore- 
grcssive process. Statistics and Probability Letters, 30:235-240, 1996. 

P.H. Franses and D. Dijk. The forecasting performance of various models for seasonality 
and nonlinearity for quarterly industrial production. International Journal of Forecast- 
ing, 21:87-102, 2005. 

P.H. Franses and M. Ooms. A periodic long-memory model for quarterly UK inflation. 

International Journal of Forecasting, 13:117-126, 1997. 
W. A. Gardner. Measurement of spectral correlation. Signal Processing, 34 (5):1111-1123, 

1986. 

W. A. Gardner, A. Napolitano, and L. Paura. Cyclostationarity: Half a century of 

research. Signal Processing, 86:639-697, 2006. 

E. G. Gladyshev. Periodically correlated random sequance. Sov. Math., 2:385-388, 1961. 
V. Gomez. The use of butterworth filters for trend and cycle estimation in economic time 

series. Journal of Business and Economic Statistics, 19(3):365-373, 2001. 
V. Gomez. Three equivalent methods for filtering finite nonstationary time series. Journal 

of Business and Economic Statistics, 17(1):109-117, 1999. 
M. Gradzewicz, J. Growiec, J. Hagemcjcr, and P. Popowski. Cykl koniunkturalny w polsce 

- wnioski z analizy spektralnej. Bank i Kredyt, 41:41-76, 2010. 
J.D. Hamilton. A new approach to the economic analysis of nonstationary time series 

and the business cycle. Econometrica, 57(2):357-384, 1989. 
A.C. Harvey and A. Jaeger. Detrending, stylized facts and the business cycle. Journal of 

Applied Econometrics, 8:231-247, 1993. 
J.R. Hicks. A Contribution to the Theory of the Trade Cycle. Clarendon Press, Oxford, 

1950. 

R.J. Hodrick and E.G. Prescott. Postwar U.S. Business Cycles: An Empirical Investiga- 
tion. Journal of Money, Credit and Banking, 29:1-16, 1997. 

H. Hurd. Correlation theory of almost periodically correlated processes. J. Multivariate 
Anal., 37:24-45, 1991. 

H. Hurd. Representation of strongly harmonizable periodically correlated process and 

their covariances. J. Multivariate Anal, 29:53-67, 1989. 
H.L. Hurd and A.G. Miamcc. Periodically Correlated Random Sequences: Spectral Theory 

and Practice. Wiley, Hoboken, New Jersey, 2007. 
R.G. King and S.T. Rebelo. Resusticating real business cycles. In: J.B. Taylor, M. 

Woodford (Eds.), Handbook of Macroeconomics, 1999. 
H-M. Krolzig. Markov-Switching Vector Autoregressions, Modelling, Statistical Inference 

and Application to Business Cycle Analysis. Springer- Verlag, Berlin, 1997. 



18 



L. Lenart. Asymptotic distributions and subsampling in spectral analysis for almost 
periodically correlated time series. Bernoulli, 17(1):290-319, 2011a. 

L. Lenart. Non-parametric frequency identification and estimation in mean function for 
almost periodically correlated time series, in preparation, 2011b. 

L. Lenart, J. Leskow, and R. Synowiecki. Subsampling in testing autocovariance for 
periodically correlated time series. J. Time Ser. Anal, 29(6), 2008. 

R.E. Lucas. Understanding business cycles. In Brunner, Karl and Meltzer, Allan H. 
(eds.), Stabilization of the Domestic and International Economy, Carnegie- Rochester 
Conference Series on Public Policy, no. 5. Reprinted in Kydland, Finn E. ( ed.), Busi- 
ness Cycle Theory, Aldershot: Edward Elgar: pp. 85-107, pages 85-107, 1977. 

S. Makridakis, S.C. Wheelwright, and R.J. Hyndman. Forecasting: Methods and Appli- 
cations. Wiley, New York, 1998. 

A. Maravall and Ana del Rio. Time aggregation and the Hodrick-Prescott filter. Banco 
de Espana Servicio de Estudios Documento de Trabajo n." 0108, 2001. 

C. Milas, P.A. Rothman, and D. Dijk. Nonlinear Time Series Analysis of Business Cycles. 
Elsevier, Amsterdam, 2006. 

W.C. Mitchell. Business cycles as revealed by business annals. In: W.L. Thorp (Ed.), 

Business Annals, 1926. 
A. Napolitano and CM. Spooner. Cyclic spectral analysis of continuous-phase modulated 

signals. IEEE Trans. Signal Process., 49 (l):30-44, 2001. 

D. R. Osborn and J. P. Smith. The performance of periodic autoregressive models in 
forecasting seasonal U.K. consumption. J Bus. Econ. Stat., 9:117-127, 1989. 

E. Parzen and M. Pagano. An approach to modeling sezonally stationary time-series. J. 
Econometrics, 9:137-153, 1979. 

D. Politis, J. Romano, and M. Wolf. Subsampling. Springer- Verlag, New York, 1999. 
M. B. Priestley. Spectral Analysis and Time Series. Academic Press, London, 1981. 

E. Serpedin, F. Panduru, L Sari, and G. B. Giannakis. Bibhography on cyclostationarity. 
Signal Processing, 85:2233 - 2303, 2005. 

P. Skrzypczynski. Analiza synchronizacji cykli koniunkturalnych w strefie euro. Materialy 

i Studia Narodowego Banku Polskiego, Zeszyt nr 210, 2006. 
P. Skrzypczynski. Wahania aktywnoci gospodarczej w Polsce i strefie euro. Materialy i 

Studia Narodowego Banku Polskiego, Zeszyt nr 227, 2008. 
J.H. Stock and M.W. Watson. Business cycle fluctuations in us macroeconomic time 

series. In: J.B. Taylor, M. Woodford (Eds.), Handbook of Macroeconomics, 1999. 
A. M. Walker. On the estimation of a harmonic component in a time series with stationary 

independent residuals. Biometrika, 58(l):21-36, 1971. 



19 




'95 '96 97 '98 '99 '00 01 '02 '03 '04 05 '06 '07 '08 '09 '95 '96 '97 '98 '99 '00 '01 '02 '03 '04 '05 '06 '07 08 '09 



(a) (b) 

Fig. 1: (a) Industrial production index in Poland (2005 year = 100%) from January 1995 to 
December 2009; (b) Realization of centered moving average filter 2xl2MA applied for industrial 
production index in Poland. 




'95 '96 97 '98 '99 '00 '01 '02 03 04 05 '06 07 '08 09 '95 '96 '97 '98 '99 00 01 02 03 04 '05 06 07 08 09 



(a) (b) 

Fig. 2: (a) First difference of centered moving average filter 2 x 12MA applied for industrial 
production index; (b) First difference of centered moving average filter 2 x 12MA applied for 
logarithm of industrial production index. 




0.44 0.87 1.31 1.75 2.18 2.62 3.05 



Fig. 3: Frequency identification using statistics n^d-i/^}) = A/n|f„(^)| and corresponding critical 
value Cn,fe(0.99%) for the realization of time series {Xt : t £ Z}: continuous line - the value of 
test statistics n„({^}) = ^/n\fn{tp)\ for ip from the set {(/c — l)7r/720 : k = 1,2, . . . , 720}; dashed 
line - critical value 41,^(99%) for V from the set {{k - l)7r/720 : A; = 1, 2, ... , 720}. 
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(b) 



Fig. 4: Frequency identification (in the set ^p,i) and corresponding amplitude estimation: (a) 

the value of test statistics Iln{{i^}) = \/n\fn{ilj)\ (continuous hne) and critical value (a) 
(dashed line) for a G {92%, 95%, 99%} and V from the set {(/c - l)7r/720 : k = 1, 2, . . . , 100}; (b) 
estimated amplitude corresponding to estimated frequencies from the set ^'p^i: X - estimated 
length of the cycle, Y - estimated amplitude. 



The value of frequency 
estimator 


i>n,l = 0.062 


^n,2 = 0.153 


'4^n,3 = 0.258 


Corresponding length 
of the cycle 
(in years) 


8.5 


3.4 


2 



Tab. 1: Estimated frequencies with corresponding length of the cycle for industry production 
index in Poland. 

A=5500 
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A=32000 
A=55000 
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Fig. 5: Business cycle in industrial production after logarithm and application of centered 
moving average filter extracted by HP filter for A = 5 500 (continuous line) A = 12 000 (dotted 
line) A = 32 000 (doted and dashed line) A = 55 000 (dashed hne). 
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Expansion 


... - Dec. 97 


Feb. 99 - May 00 


Sept. 02 - Mar. 04 


Jun. 05 - Jan. 08 


Apr. 09 - ... 


Recession 


Dec. 97 - Feb. 99 


May 00 - Sept. 02 


Mar. 04 - Jun. 05 


Jan. 08 - Apr. 09 





Tab. 2: Expansions and recessions in industrial production index in Poland in the period July 
1995 - June 2009. 



B-D_F - Mining and quarrying; manufacturing; eiectricity, gas, steam and air conditioning supply; construction 

MIGJNG.CAG - MIC Intermediate and capital goods 

MIG.ING - MIG - MIG - Intermediate goods 

MIG.CAG - Capital goods 

MIG.DCOG - MIG - Durable consumer goods 

MIG_NDCOG - MIG - Non-durable consumer goods 

B - Mining and quarrying 

C - Manufacturing 

C10-C12 - Manufacture of food products; beverages and tobacco products 

C10_C11 - Manufacture of food products and beverages 

CIO - Manufacture of food products 

Cll - Manufacture of beverages 

C12 - Manufacture of tobacco products 

C13_C14 - Manufacture of textiles and wearing apparel 

C15 - Manufacture of leather and related products 

C16 - Manufacture of wood and of products of wood and cork, except 

furniture; manufacture of articles of straw and plaiting materials 

C17 - Manufacture of paper and paper products 

CIS - Printing and reproduction of recorded media 

C19 - Manufacture of coke and refined petroleum products 

C20 - Manufacture of chemicals and chemical products 

C21 - Manufacture of basic pharmaceutical products and pharmaceutical preparations 

C22 - Manufacture of rubber and plastic products 

C23 - Manufacture of other non-metallic mineral products 

C24 - Manufacture of basic metals 

C25 - Manufacture of fabricated metal products, except machinery and equipment 

C26 - Manufacture of computer, electronic and optical products 

C27 - Manufacture of electrical equipment 

C28 - Manufacture of machinery and equipment n.e.c. 

C29 - Manufacture of motor vehicles, trailers and semi-trailers 

C29_C30 - Manufacture of motor vehicles, trailers, semi-trailers and of other transport equipment 
C31 - Manufacture of furniture; other manufacturing 
D - Electricity, gas, steam and air conditioning supply 



Tab. 3: Categorised indices describing changes in economic activity in sectors and subsectors 
of industrial production in Poland 
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Fig. 6: Logarithm of industrial production indices in Poland (2005 rok = 100%) in sectors and 
subsectors from January 1995 to February 2010. 
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Fig. 9: Frequency identification (in the set ^p,i): continuous line - the value of test statistics 
ILniii^}) = dashed line - critical value (o() for a G {92%, 95%, 99%} and from 

the set {(A; - l)7r/720 : A; = 1, 2, . . . , 100}. 
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10: Estimated amplitude and estimated length of the cycles connected with identified 
lencies in the set ^p,i: X - estimated length of the cycle, Y - estimated amplitude. 
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Fig. 11: Length of the cycles and amphtudes comparison: X - the index number, Y - estimated length of the cycle, Z - estimated 
amplitude. 
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Fig. 12: Business cycle (extracted by HP filter) in sectors and subsectors of industrial production 
after logarithm and application of centered moving average filter for A = 5 500 (continuous line) 
A = 12 000 (dotted line) A = 32 000 (doted and dashed line) A = 55 000 (dashed line). 
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